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Abstract 

The use of transgenic livestock is providing new metliods for obtaining pharmaceutically useful proteins. However, tlie 
protein expression profiles of the transgenic animals, including expression of nnilk fat globule membrane (MFGM) proteins, 
have not been well characterized. In this study, we compared the MFGM protein expression profile of the colostrum and 
mature milk from three lines of transgenic cloned (TC) cattle, i.e., expressing recombinant human a-lactalbumin (TC-LA), 
lactoferrin (TC-LF) or lysozyme (TC-LZ) in the mammary gland, with those from cloned non-transgenic (C) and 
conventionally bred normal animals (N). We identified 1, 225 proteins in milk MFGM, 166 of which were specifically 
expressed only in the TC-LA group, 265 only in the TC-LF group, and 184 only in the TC-LZ group. There were 43 proteins 
expressed only in the transgenic cloned animals, but the concentrations of these proteins were below the detection limit of 
silver staining. Functional analysis also showed that the 43 proteins had no obvious influence on the bovine mammary 
gland. Quantitative comparison revealed that MFGM proteins were up- or down-regulated more than twofold in the TC and 
C groups compared to N group: 126 in colostrum and 77 in mature milk of the TC-LA group; 157 in colostrum and 222 in 
mature milk of the TC-LF group; 49 in colostrum and 98 in mature milk of the TC-LZ group; 98 in colostrum and 132 in 
mature milk in the C group. These up- and down-regulated proteins in the transgenic animals were not associated with a 
particular biological function or pathway, which appears that expression of certain exogenous proteins has no general 
deleterious effects on the cattle mammary gland. 
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Introduction 

The technology of using genetically modified bovine lines to 
produce recombinant proteins in milk has flourished since the 
1 990s, based on the historic breakthrough of somatic cell cloning 
technology [1,2]. Milk has been proven to be an excellent vehicle 
for producing and delivering recombinant human proteins [3]. 
The cow mammary gland can be modified to synthesize foreign 
proteins in transgenic cloned bovine milk, which may be useful for 
human consumption and assimilation because of its similarities of 
composition to human milk [4]. The advantages of mUk-based 
system also include the easiness of animal housing and mainte- 
nance and lower cost of harvesting proteins than the cell 
bioreactor [5]. However, we have littie information about how 
transgenic modification and cloning process influence bovine 
endogenous gene expression, milk characteristics and health. A 
fiarther consideration is how these changes of products from these 
animals might affect the health of consumers. To address these 
issues, we examined both the protein expression profiles and 



composition of milk from transgenic cloned cows compared with 
conventionally bred animals. Many works have been done to 
analyze the compositions of meat and milk from transgenic cloned 
cattle, but the majority of these studies have focused on a few 
major constituents rather than a comprehensive analysis of 
proteins expression, especially of low-abundance proteins [6-8] . 

The concerns about consumption of food from transgenic 
cloned animals were the hot spot in recent research. The FDA has 
raised the instructions for the food from cloned animal and the 
studies have showed that the milk and meats from cloned cattle 
were as safe as the conventionally bred cattle. However, there were 
no instructions for food products from transgenic cloned animals, 
we have to use the instructions of cloned animals to examine the 
food products from transgenic cloned cattle. Milk proteins can be 
divided into three classes, namely casein, whey proteins, and 
MFGM proteins. Compared to other classes, MFGM proteins are 
the least abundant, making up 1-4% of total milk proteins, but 
they include thousands of different proteins [9], which were 
important composition of the milk proteins. 
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The MFGM is a protein-rich lipid bilayer that surrounds lipids 
in milk. Although the MFGM has been studied for 50-60 years, 
previous studies have focused primarily on membrane globule 
formation, intracellular transition, and secretion [10]. It is believed 
that the MFGM originates from the apical plasma membrane of 
mammary epithelial cells [11,12]. Precursor microlipid droplets 
from the endoplasmic reticulum fuse to each other and travel to 
the apical cytoplasm where they are surrounded by the apical 
plasma membrane and then are secreted into the alveolar lumen 
[10,13-15]. Therefore, we can obtain information about mam- 
mary epithelial cell health by analyzing the protein composition of 
the MFGM. 

In addition, MFGM proteins have many beneficial bioactivities, 
e.g., as antibiotics and anticancer agents, MUCl is one of the most 
abundant proteins in the MFGM and is thought to protect 

exposed surfaces from physical damage and invasive pathogenic 
microorganisms [16,17]. For example, bovine lactophorin C- 
terminal fragment and PAS6/7, both MFGM constituents, inhibit 
replication of human rotavirus and prevent gastroenteritis [18]. 
MFGM proteins also have commercial value in cheese and butter 
production. 

Technological advances in proteomics, especially quantitative 
proteomics, have enabled an increased understanding and 
characterization of milk: proteins [19]. So far, proteomics has 
been successfully applied to studies of human milk whey proteins, 
bovine milk whey proteins, human MFGM proteins, bovine 
MFGM proteins, and goat MFGM proteins [20-23]. 

Wc gcnc'ticalK' engineered three transgenic Ixn'ine lines that 
specifically express recombinant human a-lactalbumin (TC-LA 
group), lactoferrin (TC-LF group), and lysozyme (TC-LZ group) 
[24—26], from which the compositions of the milk include the 
human proteins with higher expression level compared to that 
from conventionally bred animals. Prior research for milk whey 
proteins and the nutritional components from transgenic cloned 
catde have indicated that the c-xpression of exogenous proteins did 
not significandy change the milk whey protein profile, and the 
mean values for the majority of the measured parameters were 
within the normal range [27]. 

In this work, we compared the MFGM proteins of colostrum 
and mature milk from these three healthy transgenic cloned lines 
to those from cloned and conventionally bred animals. The 
MFGM proteins of the five groups were characterized using 2D- 
PAGE and 2D LC-MS/MS coupled with the iTRAQ, proteomic 
strategy. Our objectives were to: (1) characterize the proteins of the 
MFGM in transgenic cloned catde; (2) compare tlu- relative 
expression of MFGM proteins among transgenic c:loned catde 
milk; and (3) investigate the health of mammary epithelial cells in 
the transgenic cloned catde. The study will supplement the former 
study by c:ompleting the evaluation for the milk products from 
transgenic cloned cattle and revealing the situation of the 
mammar)' epithelial cells from transgenic cloned catde. This 
study provides the complement evidence of the mUk composition 
analysis for the transgenic cloned bovine, and the results may 
indirectiy reflect the health of the mammary epithelial cells of the 
transgenic cloned bovine. 

Materials and Methods 

Milk Sample Collection 

The protocol was approved by the Institutional Animal Care 
and Use Committee of China Agricultural University (ID: 
SKLAB-2010-05-01). The three transgenic cattle lines TC-LZ 
(n= 10), TC-LA (n = 4), and TC-LF (n = 3) have been described 
[24—26] and were compared to cloned cattie (C, n = 3) which were 



also cloned by somatic cell nuclear transfer (SCNT), and 
conventionally bred catde (N, n = 9) of a similar genetic 
background. The cattie were similar in age and lactation period 
and were housed under the same conditions. The milk collection 
were carried out as described previously [27]. The colostrum was 
obtained during the initial three days of lactation and mature mUk 
was obtained on the 30*, 60* and 90* day after lactation and the 
milk was collected from the cows twice daily and was pooled to 
form one daily sample. Milk samples were collected by milking 
machine, and samples were immediately stored at — 20°C until 
further analysis. 

Protein Extraction 

MUk sample (50 ml) from each cow was centrifuged at 4000 xg 
to separate mUk fat globule cream from whole mUk. The separated 
cream was rinsed in 25 ml phosphate buffered saline twice and 
25 ml deionized water once. Then the milk fat globule proteins 
were obtained by methanol/trichloromethane precipitation as 
follows. The cream (2 ml) was then stored at -80°C, and upon 
thawing 10 ml methanol and 10 ml trichloromethane were added. 
The samples were centrifuged at 50,000 xg for 30 min at 4°C 
(Avanti J-26XP, Beckman Coulter, Indianapolis, USA), and the 
supernatant was discarded. The precipitated proteins were 
solubUized in a solution containing 7 M urea, 2 M thiourea, 
0.025 M DTT, 2% (w/v), and 1% CHAPS [20]. The relative 
concentration of protein in each extract was determined by 
measuring the amount of whey protein with the 2D-Quant 
quantitate kit (GE Healthcare). MFGM protein samples from each 
cow were pooled together to generate one individual sample and 
then individual samples in the same groups (TC-LZ, TC-LA, TC- 
LF, C and N) were pooled again to generate one group sample 
according to the equal protein mass. 

In-solution digestion 

Soluliic proteins (200 (Xg) were digested in-solution using the 
filter aided sample preparation (FASP) method [28] . Each protein 
extract (30 [ll) was mixed with 200 |xl of buffer containing 8 M 
urea in 0.1 M Tris-HCl (pH 8.5) in the filter unit and centrifuged 
at 1 4,000 Xg for 40 min and repeated one timt-. The flow-through 
was discarded, and iodoacetamide solution was added to the filter 
and mixed at 600 rpm for 1 min then incubated without mixing 
for 5 min. The filtered units were centrifuged at 14,000 xg for 
30 min, then 100 (tl of buffer containing 8 M urea in 0.1 M Tris- 
HCl (pH 8.0) was added to the filtered unit, then the solution was 
centrifuged again at 1 4,000 xg for 40 min and repeated twice. We 
then added 40 |xl of the above-mentioned buffer containing 
endoproteinase Lys-C (enzyme to protein ratio 1:50, w/w) and 
mixed at 600 rpm in thermo-mixer for 1 min. The units were 
incubated overnight then transferred to new collection tubes, and 
then 120 Ml of 0.05 M NH4HCO3 in water with trypsin (enzyme 
to protein ratio 1:100 w/w) was added and mixed at 600 rpm for 
1 min. The units were incubated at room temperature for 4 h and 
then centrifuged at 14,000 xg for 40 min. Finally, 50 jjl of 0.5 M 
NaCl was added and the units were centrifuged at 1 4,000 xg for 
20 min. 

The extracted peptides were desalted using a 1.3-ml C18 solid- 
phase extraction column (Sep-Pak® cartridge; Waters Corpora- 
tion, MUford, USA). The peptides were dried using a vacuum 
centrifuge and then resuspended in 5 mM ammonium formate 
containing 5% acetonitrUe, pH 3.0. 

LC-MS/MS Analysis 

LC-MS/MS was performed on a nano Acquity UPLC system 
(Waters Corporation) connected to an LTQ_ Orbitrap XL mass 
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spectrometer (Thermo Electron Corp., Bremen, Germany) 
equipped with an online nano-electrospray ion source (Michrom 
Bioresources, Auburn, USA). The MS was operated in data- 
dependent mode to switch automatically between MS and MS/ 
MS acquisition. Survey full-scan MS spectra (m/ z 300-1800) were 
acquired in the Obitrap with one microscan and with a mass 
resolution of 60,000 at m/z 400, followed by MS/MS of the eight 
most-intense peptide ions in the LTQ analyzer. LC-MS/MS 
spectra were acquired using SEQUEST (v.28 revision 12, Thermo 
Electron Corp) against the International Protein Index bovine 
database version 3. ,54. A decoy database containing the reverse 
sequences was appended to the database. The search parameters 
were set as follows: partial trypsin cleavage with two missed 
cleavage sites was considered; the peptide mass tolerance was 
1.4 Da; and the fragment ion tolerance was 1 Da. Trans 
Proteomic Pipeline software (revision 4.20; Institute of Systems 
Biology, Seattie, WA) was then used to identify proteins based on 
corresponding peptide sequences with ^95"/n confidence. The 
peptide results were filtered by PeptideProphet with a p value > 
0.90 and a ProteinProphet probability = 0.95 [29,30]. 

Quantitative Analysis of MFGM Proteins by ilRAQ LC-MS/ 
MS 

The digested peptides of the fi\'e groups were transferred to vials 
containing individual iTRAQ rc'ag(;iits (Appliixl Biosystcms) 
following the iTRAQ^ standard protocol for the 8-plex kit. The 
N group was labeled with iTRAC)113, die TC-LF with 
iTRAQ,114, TC-LZ with iTRAQ^lIS, TC-LA widi iTRAQ.116, 
and the C group witii iTRAQ_117. 

The iTRAQjlabeled samples were pooled and the SCX HPLC 
experiment was performed on a Shimadzu 20AD ,5-|j.m SCX 
polysulfoethyl column (2.1 mm x 100 mm, The Nest Group, Inc., 
MA) as the first dimension. Each coUci tc-d components of the 
processed SCX fractions ran via RP-LC ESI-MS/MS on an 
Applied Biosystems Q_-Star Elite XL mass spectrometer in which 
die RPLC column was a ZORBAX 300SB-C18 (5 nm, 
0.1 mm X 150 mm, Microm, USA). The Q-Star Elite XL mass 
spectrometer was operated in the smart information-dependent 
acquisition activated mode with automatic collision energy and 
automatic MS /MS accumulation. Survey fLdl-scan MS spectra 
(m/z 400-1800) were acquired with one microscan and a mass 
resolution of 60,000 at m/z 400, followed by MS/MS of the four 
most-intense peptide ions in the analyzer. The relative abundance 
of the MFGM proteins in the different samples was derived from 
the ionic peak areas of the iTRAQ^ reporter. 

iTRAQ_ identification and quantification analysis of the MFGM 
proteins were obtained using Protein Pilot 3.0 (Applied Biosys- 
tems, USA) with the following user-defined parameters: Sample 
Type, iTRAQ, 8-plex (Peptide Labeled); Cysteine alkylation, 
MMTS; Digestion, Trypsin; Specify Processing, Quantitate; 
Database, IPI v3.62 bovine; Search Effort, thorough ID; Results 
Quahty, Detected Protein Threshold [Unused ProtScore (Conf)] 
>1.3 (95%); Run False Discovery Rate Analysis. 

Results 

Overview of MFGM Proteome in the TC Groups 

MFGM proteins were characterized in three TC bovine lines 
and in C group and N group using 2D nano-LC-MS/MS and a 
quantitative proteomics method (iTRAQ). LC-MS/MS identified 
1225 proteins among the five groups of MFGM samples, 939 
proteins in colostrum, and 910 proteins in mature milk (Table 1). 
There were 637 proteins that were present in TC-LA group, 721 
in the TC-LF group, 720 in the TC-LZ group, 527 in the C group, 



and 668 in the N group. iTRAQ) identified 851 MFGM proteins in 
colostrum and 775 MFGM proteins in mature milk of all the 
bovine lines. The identified proteins from the three TC groups 
were compared to those from the C and N groups. The TC-LA 
group contained 166 proteins that were not found in the control 
lines, the TC-LF group contained 265 pr()t(;ins, and the TC-LZ 
group contained 184 proteins (Figure 1). We found that 43 of these 
proteins were shared by all the three TC groups (Table 2). The 
molecular mass of this subset was between 10 kDa and 90 kDa, 
and the pi values were between 4 and 10. DAVID Bioinformatics 
Resources was used to analyze the function of the 43 common 
proteins in GO terms [31]. We found that these proteins were 
enriched in the biological processes of protein transport, cellular 
component of the Golgi apparatus, and the molecular function of 
GTP binding. However, the functions of these specific expressed 
proteins did not cluster in any particular pathway (Figure 2). 

Quantitative Comparison of MFGM Proteins 

We obtained relative MFGM protein quantity information from 
the five groups using iTRAQ. Only proteins with consistent 
quantitative results in two independent iTRAQ^runs were used for 
further analysis of expression and function. A total of 459 proteins 
from colostrum MFGM and 426 proteins from mature MFGM 
were identified. Statistical analysis of expression levels in the three 
TC groups and C group compared with the N group revealed 
more than twofold differences (p^O.05) in the expression of many 
of these proteins (Figure 3). In this regard, 126 proteins in the 
colostrum and 77 proteins in the mature milk of the TC-LA group 
had a greater than twofold difference in expression, likewise, in the 
TC-LF group there were 157 proteins in colostrum and 222 
proteins in mature mUk, in the TC-LZ group there were 49 
proteins in colostrum and 98 proteins in mature milk, and in the 
cloned control group there were 98 proteins in colostrum and 132 
proteins in mature milk. 

Comparison l)et\ve(;n the MFGM proteins among the five 
groups also revealed many difierences. The unsupervised hierar- 
chical cluster analysis (HCA) and the principal component analysis 
(PCA) were used to assess the variability among the 5 groups. We 
found that the TC-LA, TC-LZ, and C groups had similar protein 
compositions in colostrum and mature milk. The TC-LF group 
was distinguishable from the other two TC groups (Figure 4). The 
colostrum of the TC-LF and TC-LA groups clustered together as 
that of the TC-LZ and C groups did. The results indicated that in 
the colostrum the TC-LZ was closer to the C group and in the 
mature milk the TC-LZ and TC-LA were closer to the C group. 
Both in colostrum and mature milk, the TC-LF groups were 
showed a greater difference from the other TC groups. In 
colostrum, 1 1 proteins were difiFerentiaUy expressed in the TC-LA 
group compared to the control groups, 25 in the TC-LF group, 
and 6 in the TC-LZ group. In mature milk, there were 8 proteins 
that were differentially expressed compared to the controls in the 
TC-LA group, 25 in the TC-LF group, and 5 in the TC-LZ group. 
In TC groups, the proteins with expression changed simulta- 
neously were not appeared. Several of these differentially 
expressed proteins, such as fibrinogen a chain, a-S2-casein, k- 
casein, P-lactoglobulin, and lactoferrin have been explored in 
other studies of whey and MFGM [20,32-34]. 

Our identified MFGM proteins were subjected to cluster 
analysis to evaluate the relationship between TC animals and 
conventionally bred controls. In the TC-LA group, six proteins 
were notably down-regulated, including a-2-antiplasmin (—4.9 
fold), fi-lactoglobulin (—4.8 fold), serum albumin (—4.9 fold), 
platelet glycoprotein 4 (—4.1 fold), lactoferrin (—4.2 fold), serine 
hydroxymethyltransferase, mitochondrial (—4.1 fold), and a 
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TC-LA(637) C+N(767) 



TC-LZ(720) C+N(767) 





TC-LF(721) C+N(767) 
D TC-LA(166) 




TC-LF(256) 



TC-LZ(184) 



Figure 1. Comparison of MFGIVI proteins from the different TC and control lines. Cloned and conventionally bred control animals [(C+N)] 
were compared with (A) MFGIVI proteins from the TC-LA line, (B) MFGM proteins from the TC-LF line, and (C) MFGM proteins from the TC-LZ line. (D) 
The intersection of differentially expressed MFGM proteins from the three TC groups yielded 43 proteins in common. 
doi:1 0.1 371/journal.pone.01 05378.g001 



Table 1 . Proteins identified in the MFGM of colostrum and mature milk of the TC-LZ, TC-LA, and TC-LF transgenic cloned lines and 
of the cloned and conventionally bred normal controls. 







TC-LZ 


TC-LA 


TC-LF 


Cloned 


Normal 


Total 


Colostrum 


620 


516 


411 


220 


562 


939 


Mature milk 


444 


393 


567 


464 


353 


910 


in total 


720 


637 


721 


527 


668 


1225 



doi:l 0.1 371 /journal.pone.Ol 05378.t001 
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Table 2. Proteins specifically expressed in transgenic cloned cattle MFGM. 




No 


Accsssion 


Description 


MW 


PeptideLencfth 


Isoelectric 
Point 


1 


IPI00686092.1 


Peroxiredoxin-1 


22194.98 


199 


8.6084 


2 


IPI00687416.2 


Rras2 Protein 


23384.41 


204 


5.8076 


3 


IPI00687550.1 


3'(2'),5'-Bisplnosphate 
Nucleotidase 1 


33306.63 


308 


5.1936 


4 


IPI00689149.2 


Oxysterol-Binding Protein 


84236.69 


746 


6.2374 


5 


IPI00692733.2 


Transmembrane Protein 30A 


40656.91 


361 


8.6073 


6 


IPI00693597.3 


Similar To Adiponutrin 


50423.45 


455 


7.5812 


7 


IPI00695308.5 


Mucin-16 (Fragment) 


31335.98 


282 


9.0396 


8 


IPI00695670.1 


Upf0585 Protein C160rfl3 Homolog 


22699.17 


204 


7.9092 


9 


IPI00695741.2 


62 Kda Protein 


62274.23 


569 


5.3891 


10 


IPI00696647.1 


Similar To Pincher Isoform 1 


60985.31 


540 


7.671 1 


1 1 


IPI00697565.3 


Ras-Related Protein Rab-21 


241 30.66 


222 


8.0481 


12 


IPI00698430.1 


Eul<aryotic Translation Initiation 
Factor 5 Isoform 7 


48939.49 


429 


5.2609 


13 


IPI00700391.1 


Similar To V-Crk Sarcoma 
Virus CtIO Oncogene 
Homolog Isoform 1 


33809.31 


304 


5.2451 


14 


IPI00703351.2 


Similar To Aldehyde Dehydrogenase 
3B1 Isoform 1 


53661.38 


486 


8.4692 


15 


IPI00703423.1 


B-Cell Receptor-Associated Protein 31 


27883.64 


245 


9.765 


16 


IPI00704666.1 


Trna Methyltransferase 112 Homolog 


14309.9 


125 


5.8553 


17 


IPI00704752.1 


Ras-Related Protein Rab-7A 


23528.72 


207 


6.5881 


18 


IPI00706451.3 


Adp-Ribosylation Factor 4 


20514.36 


180 


5.9819 


19 


IPI00706632.1 


Elongation Factor 1 -Gamma 


50345.45 


440 


6.5236 


20 


IPI00706968.3 


Nuclear Factor Of Kappa 

Light Polypeptide Gene Enhancer 

In B-Cells Inhibitor, Beta 


37543.32 


357 


4.2979 


21 


IPI00707103.4 


Synaptosomal-Associated Protein 29 


28471.85 


258 


5.1001 


22 


IPI00707616.3 


Cidea Protein 


24524.03 


219 


9.7273 


23 


IPI00708591.1 


Znrf2 Protein 


23870.75 


238 


7.1119 


24 


IPI0071 0727.1 


Transitional Endoplasmic Reticulum Atpase 


89272.24 


806 


4.8908 


25 


IPI0071 1 304.4 


Similar To Inositol 
Polyphosphate-5-Phosphatase A 


70475.79 


637 


8.6043 


26 


IPI0071 3743.1 


lnosine-5 '-Monophosphate 
Dehydrogenase 2 


55726.54 


514 


7.339 


27 


IPI00714515.1 


Fgr Protein 


59292.95 


527 


5.3694 


28 


IPI00714621.3 


Uncharacterized Protein 
C10rf93 Homolog 


21482.83 


201 


6.2609 


29 


IPI00714818.1 


Eukaryotic Translation Initiation 
Factor 3 Subunit H 


39879.94 


352 


6.653 


30 


IPI0071 4992.4 


Cgmp-Dependant Type li Protein Kinase 


87033.26 


762 


8.6701 


31 


IPI00715153.2 


Fas-Associated Factor 2 


52629.57 


445 


5.4387 


32 


IPI0071 6638.7 


Wd Repeat-Containing 
Protein 1 


66214.06 


606 


6.7147 


33 


IPI0071 7376.3 


Similar To Cytoskeleton-Associated 
Protein 4 


64489.56 


589 


5.9276 


34 


IPI00718291.4 


Rab14 Protein 


23881.49 


215 


6.1266 


35 


IPI00718671.5 


Synaptosomal-Associated Protein 


23212.11 


211 


4.4574 


36 


IPI00729755.2 


Acyl-Coa Synthetase Long-Chain 
Family Member 1 


78228.82 


699 


7.2276 


37 


IPI00730045.3 


Similar To Rififylin Isoform 4 


39655.22 


356 


6.1033 


38 


IPI00733615.1 


Dnaj Homolog Subfamily 
A Member 2 


45731.44 


412 


6.4421 


39 


IPI00867128.1 


Phosphoserine Aminotransferase 


40502.07 


370 


7.7533 



PLOS ONE I www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | e105378 



Comparative Proteomics of MFGM Proteins from TC Cattle 



Table 2. Cont. 





No 


Accession 


Description 


MW 


PeptideLength 


isoelectric 
Point 


40 


IPI00867311.1 


Scamp2 Protein 


36678.97 


328 


6.1107 


41 


IPI00906554.1 


Guanine Nucleotide-Binding 
Protein Subunit Alpha-1 1 


42042.75 


359 


6.029 


42 


IPI00924133.1 


Similar To Ciliary Neurotrophic 
Factor Receptor Alpha Precursor 
(Cntfr Alpha) Isoform 3 


40661.3 


372 


6.9436 


43 


IPI00944416.1 


CybSB Protein 


16267.5 


146 


4.8233 



doi:l 0.1 371 /journal.pone.OI 05378.t002 



putative uncharacterized protein (—4.4 fold). In the TC-LF group, 
10 proteins were notably different from the control group, 
including apolipoprotein A-II (—4.4 fold), ATP-binding cassette 
sub-famUy G (WHITE) member 2 (+4.3 fold), a- 1 -acid glycopro- 
tein (—4.9 fold) similar to solute carrier family 39 (zinc transporter) 
member 8 (+4.1 fold), ot-S2-casein (+4.2 fold), MUG 1 (+4.9 fold), 
acyl-GoA synthetase long-chain family member 1 (+4.8 fold), and 
vitamin D-binding protein (—5.2 fold). In the TG-LZ, four 
proteins were notably different from the control group, namely 
fibrinogen a chain (+2.8 fold), P-lactoglobulin (—2.7 fold), cidea 
protein (-3.1 fold), and similar to ATPase type 13A4 (-2.9 fold). 
In the C group, six proteins were notably different from the 
conventionally bred group, namely fibrinogen a chain (+3.1 fold), 
a-S2-casein (+5.3 fold), a-Sl-casein (+4 fold), Ritl protein (-3.6 
fold), ot-lactalbumin (—4.8 fold), and nucleobindin-1 (+3.3 fold). 

Function and Pathway Analysis of Differentially 
Expressed Proteins 

The Ingenuity Pathway Analysis (IPA, Ingenuity Systems, Inc. 
Redwood City) software was used to analyze the functions and 
networks of the differentially expressed proteins. To further 
examine the specifically expressed proteins, they were categorized 
according to their biological process, cellular component or 
molecular function as annotated in the Gene Ontology (GO) 
database. Differentially expressed MFGM proteins from TG-LA 
colostrum were mainly enriched in the biological processes of 
response to wounding and the acute inflammatory response, in the 
ceUular components of the extracellular region, and in the 
molecular function of endopeptidase inhibitor activity. Differen- 
tially expressed MFGM proteins from TG-LA mature milk were 



mainly enriched in cytoskeleton organization, membrane-bounded 
vesicles, and unfolded protein binding. Differentially expressed 
MFGM proteins from TG-LF colostrum were mainly enriched in 
response to organic substance, the extracellular region, and 
endopeptidase inhibitor activity. Differentially expressed MFGM 
proteins from TG-LF mature milk were mainly enriched in smaU 
GTPase-mediated signal transduction, pigment granules, and 
GTPase activity. DiflferentiaUy expressed MFGM proteins from 
TG-LZ colostrum were mainly enriched in response to wounding, 
the extracellular region, and cell-surface binding. Finally, differ- 
entially expressed MFGM proteins from TG-LZ mature milk were 
mainly enriched in cell motility, the actin cytoskeleton, and actin 
binding (Figure 5). 

In the TG-LA group, there were 18 proteins that were changed 
above 3 folds in colostrum and mature mUk. The pathway analysis 
showed that the changed MFGM proteins in the TG-LA 
colostrum most involved in the network which have the functions 
of cellular movement, immune cell trafficking, ceU-to-ceU signaling 
and interaction. And in the mature milk are cell death, ceU-to-ceU 
signaling and interaction, hematological system development and 
function (Figure 6A, B). 

In TG-LF group, there were 46 proteins that were changed in 
colostrum and mature milk significandy. The IPA software showed 
that the changed MFGM proteins in mature milk involved in the 
network of which top functions are cellular assembly and 
organization, cellular function and maintenance, dermatological 
diseases and conditions. And the MFGM proteins from colostrum 
involved in the network of which top functions are lipid 
metabolism, molecular transport, and smaU molecule biochemistry 
(Figure 6G, D). 




Figure 2. Functional analysis of MFGM proteins specifically expressed in the transgenic cloned animals. Analysis of proteins enriched in 
the transgenic cloned groups according to: (A) biological process, (B) cellular component, and (C) molecular function. Blue curves indicate the degree 
of enrichment of gene function. 
doi:1 0.1 371/journal.pone.01 05378.g002 
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Figure 3. Cluster analysis of MFGM proteins. (A) Analysis of MFGM proteins in colostrum of the three TC groups compared to cloned [(C)] and 
conventionally bred [(N)] controls. (B) Analysis of MFGM proteins in mature milk of the three transgenic groups compared to controls. 
doi:10.1371/journal.pone.0105378.g003 




Figure 4. Principal component analysis of MFGM proteins from the three transgenic groups, the cloned group, and the 
conventionally bred control group. (A) Analysis of MFGM proteins in colostrum. (B) Analysis of MFGM proteins in mature milk. 
doi:10.1371/journal.pone.0105378.g004 
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Figure 5. Molecular function analysis of differentially expressed MFGM proteins in the transgenic groups. Enrichment analysis of 
differentially expressed MFGM proteins in (A) TC-LA colostrum, (B) TC-LA mature milk, (C) TC-LF colostrum, (D) TC-LF mature milk, (E) TC-LZ colostrum, 
and (F) TC-LZ mature milk. Blue curves indicate the degree of enrichment of gene function. 
doi:10.1371/journal.pone.0105378.g005 




- A B C 




D E F 

Figure 6. Pathway analysis of differentially expressed proteins from the transgenic groups. IPA defined pathways of the differentially 
expressed MFGM proteins from: (A) TC-LA colostrum, (B) TC-LA mature milk, (C) TC-LF colostrum, (D) TC-LF mature milk, (E) TC-LZ colostrum, and (F) 
TC-LZ mature milk. 

doi:1 0.1 371/journal.pone.01 05378.g006 
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In group of TC-LZ, there were 10 proteins that were altered 
obviously in colostrum and mature mUk. The IPA software showed 
that the changed MFGM proteins in TC-LZ mature milk involved 
in the network of which the top functions are cell-to-cell signaling 
and interaction, tissue development, cell death. And in colostrum 
MFGM proteins, the top functions are inflammatory response, 
ceU-to-ceU signaling and interaction, cell signaling (Figure 6E, F). 

The MFGM proteins of C group were also compared with the 
N group, the results showed that the changed proteins in the 
colostrum involved in the functions of cell-to-cell signaling and 
interaction, tissue development, hematological system develop- 
ment and function. In mature MFGM, the changed proteins were 
in the network of functions: cellular movement, cellular assembly 
and organization, cellular function and maintenance. Differential- 
ly expressed MFGM proteins in colostrum from the C group were 
enriched in the biological processes of response to wounding and 
the acute inflammatory response, in the cellular components of the 
extracellular region, and in the molecular function of GTPase 
activity. Differentially expressed MFGM proteins in mature milk 
of the cloned control animals were enriched in the biological 
process of glucose catabolism, in the cellular component of the 
actin cytoskeleton, and in the molecular function of cytoskeletal 
protein binding (Figure 7). 




Discussion 

Proteomics is one of the most popular methods to analyze 
complex samples. Recent advances in proteomics make the 
identification of low-abundance proteins and analysis of complex 
proteins such as high-molecular-weight transmembrane proteins 
possible, such analyses are not possible with standard 2D-PAGE 
methods alone. The combination of LC-MS/MS and ITRAQ, or 
methods with the iTRAQ^ technologies are suitable for analysis of 
hydrophobic macro-molecules like plasma membrane proteins and 
are highly sensitive, fast, and reliable. Comparing with traditional 
2D-PAGE method, the methods based on mass spectrometry are 
more suitable for MFGM proteins analysis. Data from several 
previous studies of MFGM proteins in humans [23], cattle [35] 
and goats [36] demonstrated that the amount of the identified 
proteins by this method was larger than other methods and the 
iTRAQ^ method was suitable to analyze the MFGM proteins with 
the same sensitivity compared to LC-MS/MS. 

Applying these methods, we examined the effects of exogenous 
gene expression and cloning techniques on bovine MFGM 
proteins, an neglected area of research. Comparing three TC 
bovine lines expressing different exogenous proteins with C and N 
group, we identified 43 proteins that were specifically expressed in 
the transgenic animals below the detection limit for 2D-PAGE and 
silver staining. Considering both the concentration and functional 




Figure 7. Molecular function and pathway analysis of differentially expressed MFGM proteins from the cloned control group. (A, Bj: 

IVlolecular function of IVIFGIVl proteins from colostrum (A) and mature milk (B). (C, D): IPA defined pathways of MFGIVI proteins from colostrum (C) and 
mature milk (D). 

doi:1 0.1 371/journal.pone.01 05378.g007 
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Figure 8. Pathway analysis involved in the Lipid Metabolisms of the TC-LA MFGM proteins. 

doi:1 0.1 371 /journal.pone.01 05378.g008 



analysis, we believe that these specific expressed proteins did not 
adversely influence the composition of mUk or the health of these 
transgenic animals. Peroxiredoxin- 1 was hypothesized to respon- 
sible for the anti-oxidative effect of laminar share stress [37]. 3' 
(2'), 5'-bisphosphate nucleotidase 1, inosine-5' -monophosphate 
dehydrogenase 2, elongation factor 1 -gamma, sdp-ribosylation 
factor 4 and phosphoserine aminotransferase involved in the 
metabolism of the cell, the protein ras-related protein rab-2 1 and 
guanine nucleotide-binding protein subunit alpha- 1 1 involved in 
the mediated signal transduction. The function of b-cell receptor- 
associated protein 31 involves in apoptotic cleavage of cellular 
proteins, the function of ras-related protein rab-7A involves in 
protein transport, the function of acyl-coa synthetase long-chain 
family member 1 involves in the fatty acid degradation and 
metabolism, the function of Eukaryotic Translation Initiation 
Factor 3 Subunit H was in regulation of translational initiation. 

There were distinguishable differences between the three 
transgenic lines in terms of exogenous protein expression. The 
TC-LF group expressed recombinant human lactoferrin at high 
levels, i.e., 2.5-3.8 g/1 [25], the TC-LA group, expressed 
recombinant human ot-lactalbumin at intermediate levels, i.e., 
1.5 g/1 [26], and the TC-LZ group expressed recombinant human 
lysozyme at low levels, i.e., 13-28 mg/1 [24]. The expression 
profile of TC-LF group displayed a marked difference with that of 
other TC groups but the range of variation was comparable to the 
C group. Comprehensive analysis of the relative quantities of 
MFGM proteins in the transgenic and control groups led to the 
hypothesis that MFGM protein expression was influenced by the 
quantity of the exogenous proteins. In the TC-LF group. 



lactoferrin expression was higher than that of the exogenous 
proteins in the other two transgenic groups, and the number of tlie 
affected MFGM proteins was also higher. One possible mecha- 
nism for this could entail the S 1 00 proteins, which are involved in 
inflammatory and antimicrobial responses. We found low abun- 
dance of the S 1 00A8 and S 1 00A9 proteins in colostrum of the 
transgenic animals relative to the N group. This may be a 
consequence of the fact that LF, LA, and LZ (i.e., which were 
expressed exogenously in our study animals) are all involved in the 
inflammatory response. 

The MFGM proteins are important in cattle cultivation and 
quality of livestock products, especially in cheese industry. The 
MFGM proteins were the major composition in MFGM, they play 
a core role in mammary epithelial cell against bacteria, so these 
proteins and immune globulin proteins are highly expressed in 
colostrum MFGM, and the colostrum was beneficial to the 
acquired immunity in new born. 

In TC-LA group, we analyzed the functions of the proteins and 
found out that plenty of them were not defined, such as similar to 
lipocalin, oi-2-antiplasmin, Cd5L protein, putative uncharacterized 
protein, canx protein. Hemopexinand serotransferrin coming from 
blood majorly influenced by the cattle individually. Actin-related 
protein 2/3 complex subunit 3, rab35 protein, serine hydro- 
xymethyltransferase, mitochondrial, 1-lactate dehydrogenase B 
Chain are involved in ubiquitin ligase [38]. Apolipoprotein E 
plays a key role in regulating plasma levels of lipoproteins and its 
level was not associated with high-density lipoprotein fractions 
[39]. The former studies of the composition of mUk from 
transgenic cloned cattle indicated that the fatty acids of TC-LA 
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colostrum were less than that of the N group especially in 
monounsaturated fatty acids [27] . The pathway analysis found the 
down-regulated lipid metabolism (Figue 8), which presumably was 
involved in monounsaturated fatty acids synthesis. 

In TC-LF group, some proteins liave high abundance in milk 
such as K-Casein, 0(-S2-Cascin, scrotransferrin, a-lactalbumin,- 
serum albumin and lactoperoxidase, so a small amount of residual 
proteins will cause detectable change of concentrations. Plenty 
proteins with the function of lipid metabolism were differentially 
expressed. Bovine MUC 1 prevents binding of bacteria to human 
intestinal cells and has a role in preventing the binding of common 
enteropathogenic bacteria to human intestinal epithehal surfaces 
[40]. The expression of fatty acids of mature TC-LF milk was less 
than that of other four groups, especially in saturated fatty acids 
[27]. Plenty proteins which involve in lipid metabolism were 
expressed lowly in TC-LF mature MFGM, such as apolipoprotein 
A-II, which is one of the major proteins in HDL and its main 
function is to modulate the lipid binding and lecithin-cholesterol 
acyltransferase activities of HDL by promoting the dissociation of 
apo A-I from HDL [41]. ATP-binding cassette, sub-family G 
(White), member 2 have the function of maintaining lipid 
homeostasis in the mammary gland [42]. Fatty acid-binding 
protein (FABP) regulates the charmeling of fatty acids toward 
copious milk fat synthesis in bovine mammary gland [43]. The 
expression level of CIDE-A protein was regulated by insulin and/ 
or fatty acids in mammary epithelial cells, and thereby played an 
important role in lipid and energy metabolism [44]. Acyl-coa 
synthetase long-chain family member 1 (ACSLl) which involves in 
fatty acid analysis [43] was differentially expressed in lactation,. 
ACSLl, FABP and CD 36 were suggested to have the function of 
channel long chain fatty acids toward esterification into milk 
triacylglycerol [43,45]. 
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The studies showed that both the bovine and human a- 
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In summary, the study use.s proteomics methods to analyze 
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MFGM protein profiles, and the relative quantity expression of 
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